Magneto-localization in Disordered Quantum Wires 



Stefan Kettemann, Riccardo Mazzarello 
/. Institut f. Theoretische Physik, Universitat Hamburg, Germany 

The magnetic field dependent localization in a disordered quantum wire is considered nonpertur- 
batively. An increase of an averaged localization length with the magnetic field is found, saturating 
at twice its value without magnetic field. The crossover behavior is shown to be governed both in the 
weak and strong localization regime by the magnetic diffusion length Lb- This function is derived 
analytically in closed form as a function of the ratio of the mean free path I, the wire thickness W, 
and the magnetic length Is for a two-dimensional wire with specular boundary conditions, as well as 
for a parabolic wire. The applicability of the analytical formulas to resistance measurements in the 
strong localization regime is discussed. A comparison with recent experimental results is included. 
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I. INTRODUCTION 



The phase coherent movement of electrons in a disorder potential can result in strong localization due to quantum 
interference . As soon as the localization length L c becomes smaller than the size of the sample L and the phase 
coherence length, L^, the resistance increases exponentially. 

The strong localization due to quantum interference is known to depend on the global symmetry of the disordered 
electron system In disordered quantum wires, the localization length is 



(3t{?ivSDq, 



(1) 



where (3 = 1,2,4, corresponding to no magnetic field, finite magnetic field, and strong spin- orbit scattering or 
magnetic impurities, respectively. v{E) is the electronic density of states in the wire. Do — TVp/d is the classical 
diffusion constant of the electrons in the wire, with r the elastic scattering time, vf the Fermi velocity, and d the 
dimension of classical diffusion. S is the wire crossection. This result was first obtained by calculating the spatial 
decay of the density correlation function for wires with diffuse crossections and many transversal channels N 3> 1. 
It can also be obtained by calculating the transmission probability through thin, few channel wires upto a correction 
of order l/N: L c = (f3N + 2 — 0)1 jjj, where I — vpr is the mean free path, and = 1, 2, 4, as defined above. This 
correction ensures that the localization length is for a single channel, N = 1, independent of (3, L c = 21. 

Recently, the doubling of the localization length was observed in sub-micron thin wires of Si 5- doped Ga As 
structures by Khavin, Gcrshcnson and Bogdanov, who found a continously decreasing activation energy when the 
magnetic field is increased, saturating indeed at one half of its field free value || . This symmetry dependence of the 
localization properties of quantum wires allows to test our present theoretical understanding by detailed comparison 
with the experiment. The quantum wires used in the experiment have mean free paths which are smaller than or 
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comparable to their thickness. Also, in addition to the disorder in the bulk due to the random electrostatic potential 
of the donor impurities, there is an unspecified surface roughness which may influence the classical mobility of the 
wires as well as its quantum transport properties. Therefore, a more detailed analysis of the localization length as 
function of these parameters is called for, in order to be able to compare the theory with the experimental results 
quantitatively. 

In the next section we review the known weak localization corrections to the conductivity in disordered quantum 
wires and their magneto-sensitivity as function of mobility, wire thickness, and electron density J5|-3|. 

In the third section, the non-perturbative theory of localization in disordered electron systems [2] is extended, in 
order to allow the study of wires with ballistic crossections. 

In the fourth section, the magnetic phase shifting rate is introduced and identified with a correlation function of the 
magnetic vector potential, relating it to the coefficient of the time reversal symmetry breaking term in the nonlinear 
sigma model. This expression for the magnetic phase shifting rate, is calculated anayltically for arbitrary ratios of 
the mean free path I and the width of the wire W, and compared with previously derived analytical and numerical 
results HJ^] for a wire with specular boundary scattering. 

Next, it is calculated for a wire with harmonic confinement which allows to extend the analysis to stronger magnetic 
fields, when the cyclotron radius, lc is smaller than the the wire thickness W, but still larger than the elastic mean 
free path. In that regime a new enhancement mechanism for the magnetic phase shifting rate leading to a stronger 
magneto-sensitivity, is identified. 

In the fifth section, the autocorrelation function of spectral determinants (ASD) |9|Jl0(| is considered for a coher- 
ent disordered quantum wire, which shows the expected crossover from Wigner- Dyson statistics [jll), typical for a 
spectrum of extended states in phase coherent disordered metal systems [||, to Poisson statistics, corresponding to a 



spectrum of localized states [ 12 - 17J] , as the length of the wire is increased beyond a localization length Lc, as reported 
earlier O]. 

This crossover length scale to Possionian statistics is used to derive the averaged localization length of disordered 
quantum wires, and it is shown that it yields the correct symmetry dependence, Eq. (Q). A comparison with the 
result of the supersymmetric theory of the two-terminal conductance of a disordered quantum wire, is given. It is 
concluded, that the definition of an averaged localization length, by the decay of an energy level correlation function, 
can be used to consider analytically the magnetic field dependence of the localization length. Thereby, analytical 
formulas for the localization length as a function of wire width, mean free path and magnetic field are derived. 

In the sixth section, the theory of finite tempreature magnetoresistance in quantum wires is discussed. In particular, 
the variable range hopping conductivity in quantum wires is reviewed for various temperature and dimensional regimes. 
It is shown that in a wide temperature regime the resistance has an activated behaviour, and that therefore, the 
activation gap can be directly measured and related to the localization length of the electrons in the wire. 

This allows a comparison of the analytical results for the magnetic field dependence of the localization length with 
these experimental results, as done in the seventh section. 

In appendix A, the functional integral representation of the ASD by Grassmann intergals is given, and the averaging 
over disorder is performed. In appendix B the derivation of the magnetic phase shifting rate is given. In appendic C 
the representation of the matrix fields Q is given, and their Laplacian derived. 



II. WEAK LOCALIZATION 

Classically, the transport of a disordered conductor is characterized by its mobility /i = qr/m and the electron 
density n related to the the classical Drude conductivity oq — nq 2 r jm. Alternatively, it can be characterized by the 
diffusion constant D, which is in a metal related to the conductivity by the Einstein relation do = 2q 2 vD. 

When the electrons diffuse coherently, quantum interference without magnetic field results in a suppression of the 
conductivity of a quantum wire of order p3| 



Aa 



1), (2) 



where is the phase coherence time, that increases when decreasing the temperature as a power law: 

H ~ T-i (3) 

and defines the phase coherence length, which an electron diffuses coherently, = (Dt^) 1 / 2 . 

Quasi elastic electron-electron scattering can be the dominant low temperature dephasing mechanism and yields 
7 = 2/3 for a 1-d wire and 7 = 1 for a 2-d film |7p4|]. At higher temperatures the exponent crosses over to 7 = 4 due 
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to electron-phonon scattering at temperatures /cbT <C (h 2 /t6f)^d where flo is the optical Debye phonon frequency. 
This power can be smaller, due to the confinement, in quantum wires. 

The above definition of the phase coherence rate is not applicable when approaching the localized regime, and the 
phase coherence length is larger than the localization length L c . Also, there are mechanisms which may lead to a 
saturation of below T — IK, as observed in a wide range of conductors p5| , p6| ]. 

A magnetic field breaks the time reversal symmetry. Therefore, the magnetic phase accumulated in a Brownian 
motion of electrons, enters effectively as an additive contribution to the phase coherence rate, diminishing the weak 
localization corrections of the conductivity pi]] . For wires with diffusive width W > I, it varies quadratically with the 

magnetic field, 1/t${B) = l/r^ + D^SB 2 /Kd, where S is the crossection of the wire, and the constant Kd depends 
on the geometry of the wire, the direction of the magnetic field and the scattering mechanisms H. For example, for 
a 2-dimensional wire of diffusive crossection in a perpendicular magnetic field, it yields, Kd = 3. In this way, the 
conductivity increases to its classical value, when the magnetic field is turned on. 

For a wire with ballistic crossection and a magnetic field being perpendicular to its crossection, the magnetic field 
dependence of the weak localization correction to the conductivity is weakened by flux cancellation effects due to 
boundary scattering ||. If the magnetic field is so small that less than one flux quantum (fi — h/e is penetrating an 
area Wl, the effective dephasing rate 1/t^(B) is quadratically increasing as for diffusive crossections. Its slope was 
found to be by at least a factor W/l smaller, as a consequence of the flux cancellation effect of edge to edge skipping 
orbits §,§. 

When BWl 3> <fio, the effective dephasing rate I/t^B) was found by a semiclassical method, to increase only 
linearly with the magnetic field B in this regime J6|||. 

In the presence of magnetic impurities, scattering the electrons with a rate 1/t<j, there is no temperature dependence 
of the conductivity, if 1/rg ^> 

Strong spin-orbit scattering reverses the sign of the quantum correction to the conductivity |^7j . The conductivity 
is then larger than classically expected. This can be observed by increasing an external magnetic field, which destroys 
time reversal invariance and acts through an effective decoherence time 1/t^B) = l/r^ as noted above. In the case 
of moderately strong spin-orbit scattering, the conductivity decreases therefore when the magnetic field is turned on 

§• 

At low temperatures, when the dephasing rate I/r^ becomes smaller than the typical energy scale of strong lo- 
calization, the local level spacing Ac = \/{vWLc), a perturbation theory in the elastic scattering rate I/r is no 
longer ap propr iate, and a nonperturbative treatment of disorder is called for, as the scaling theory of localization does 
indicate |l9|]2Ctl. 



III. NONPERTURBATIVE THEORY OF LOCALIZATION IN DISORDERED QUANTUM WIRES 

In this section, the nonperturbative theory of disordered noninteracting electrons in quantum wires is derived 
p8| , p2| j2]|. Its action, governed by the long wave length modes corresponding to diffusion, the nonlinear sigma model 
is rederived, extending previous derivations, to allow the description of quantum wires with ballistic crossections. 

The Hamiltonian of disordered noninteracting electrons is 

H = e (p - qA) + V(x) + <rb s (x) + au S o x p, (4) 

where q is the electron charge. In the following, we will generally approximate the electronic dispersion e (p — qA) 
by (p — qA) /(2m), where m is the effective electron mass, but note that higher moments are sometimes needed to 
regularize the correlation functions, calculated below. 

V^(x) is taken to be a Gaussian distributed random function (V^(x)) = 0, and (V^(x)V^(x')) = TiASL/(2ttt)5(x — x'), 
which models randomly distributed, uncorrelated impurities in the sample. A = 1/{vSL) is the mean level spacing. 
This corresponds to a Gaussian distribution function P(V) = cxp(— J -^-V(x) 2 ) of the disorder potential, defining 
the disorder average as < ... >v— J Tlx dVP(V).... According to the central limit theorem, this is therefore a good 
description of the various sources of randomness in the electrostatic potential, in which the electrons are moving. 

The vector potential is used in the gauge A = (—By, 0,0), where x is the coordinate along the wire of length L, 
y the one in the direction perpendicular both to the wire and the magnetic field B, which is directed perpendicular 
to the wire. The angular brackets denote averaging over impurities, a is the electronic spin operator, and b s (x) is a 
random magnetic impurity field. Ugo is the local electrostatic field of impurities with large atomic number Z, which 
do give a stronger spin orbit coupling to the conduction electrons. 

The Hamiltonian can be classified by its symmetry with respect to time reversal and spin rotation as summarized 
in Table I. 
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It has been noted that the averaged density of states or the averaged one-particle Green's function does not contain 
any information on the localization of Eigenfunctions of the disordered Hamiltonian H [p8[ . The physical reason 
is, that the one-particle Green's function describes the propagation of the wave function amplitude ^(x). Elastic 
impurity scattering randomizes the phase of the amplitude and therefore, this propagator decays on the scale of the 
mean free scattering time t. To catch classical diffusion and quantum localization, at least the evolution of the density 
or amplitude square has to be averaged over the disorder, leading to a correlation function of two one-particle Green's 
functions. While weak localization corrections can be calculated within a diagrammatic perturbation expansion of 
such correalation functions [f|[23), the study of strong electron localization in a disordered potential, necessitates 
a nonperturbative averaging of such products of Green's functions. This can be achieved by means of the super- 
symmetry method, whereby the product of Green's functions is written as a functional integral §). Thus, the average 
over the form of the disorder potential can be done right at the beginning as a Gaussian integral, exactly. 

Here, for simplicity, we present the derivation of a simpler correlation function, which does not necessitate the use 
of the full super-symmetry method, but still contains some information on strong quantum localization, as shown 
recently (l8H|,p 



The statistics of discrete energy levels of a finite coherent, disordered metal particle is an efficient way to characterize 
its properties Q . This can be studied by calculating a disorder averaged autocorrelation function between two energies 
at a distance u> in the energy level spectrum. Thereby, an uncorrelated spectrum of localized states can be distinguished 
from a correlated spectrum of extended states. 

The autocorrelation function of spectral determinants (ASD) is the most simple such spectral correlation function, 
which allows to explore complex quantum systems analytically, and still does contain nontrivial information on level 
statistics and, thus, on localization Jlsjjz^l . 

It is an oscillatory function whose amplitude decays with a power law, when the energy levels in the vicinity of the 
central energy E are extended, while a Gaussian decay is a strong indication that all states are localized. 

It is defined by C{u) = C(u)/C(0), C{lo) = (det(£ + u/2 - H)det(E - u/2 - H)) , where E is a central energy. 

Since it is a product of two spectral determinants, and a spectral determinant can be written as a Gaussian functional 
integral over Grassmann variables ip, ip*, one does need at least a 2-component Grassman field, one for each spectral 
determinant. 

In general, 4a -component Grassman fields arc needed to get the functional integral representation of the ASD. Here, 
a — 1, when the Hamiltonian is independent of the spin of the electrons, and each level is doubly spin degenerate. 
There is one pair of Grassman fields for each determinant in the ASD and each pair is composed of a Grassman field 
and its time reversed one, as obtained by complex conjugation, a — 2 has to be considered, when the Hamiltonian does 
depend on spin, as for the case with moderately strong magnetic impurity or spin- orbit scattering. This necessitates 
the use of a vector of a spinor and the corresponding time reversed one. 

The representation as a Gaussian functional integral over Grassmann variables is given explicitly for a = 1 in 
appendix A. There, the averaging over disorder and the decoupling of the resulting ip 4 interaction with a Gaussian 
integral over a matrix field Q is given. Thus, the disorder averaged ASD is given by a functional integral over a matrix 
field Q. 

The matrix Q is element of the full symmetric space, including rotations between the subspace corresponding to the 
left and the right spectral determinant. Therefore, the long wavelength modes of Q, do contain the nonperturbative 
information on the diffuson and Cooperon modes. 

In order to consider the action of long wavelength modes governing the physics of diffusion and localization, one 
can now expand around the saddle point solution of the action, satisfying for u) = 0, 

Q = < x | 1/(E -H + iT%/{2t)Q) I x > . (5) 

This saddle point equation is found to be solved by Qo = A. For a = 1, and B = 0, at ui = 0, the rotations U, 
which leave Q in the symplectic symmetric space yield the complete manifold of saddle point solutions as Q — UAU, 
where UU = 1, with Q T C — CQ. The modes which leave A invariant, elements of Sp(l) x Sp(l) are surplus, or 
spontanously broken, and can be factorized out, leaving the saddle point solutions to be elements of the symmetric 
space Sp(2)/(S P (1) x Sp(l)) 0. 

For a — 2 the matrix C is, due to the time reversal of the spinor, substituted by ia%r\ p^ |. Both magnetic 
impurities and spin-orbit scattering reduce the Q matrix to unity in spin space. Thus, C has effectively the form t\. 
The condition Q T C — CQ leads therefore to a new symmetry class, when the spin symmetry is broken but the time 
reversal symmetry remains intact. This is the case for moderately strong spin-orbit scattering. Then, Q are 4x4- 
matrices on the orthogonal symmetric space 0(4)/(0(2) x 0(2)) ^8|, which is the nonperturbative consequence of 
the sign change of a spinor component under time reversal operation, which leads to the positive quantum correction 
to the conductivity in perturbation theory p3] . With magnetic impurities both the spin and time reversal symmetry 
is broken, and the Q- matrices are in the unitary symmetric space U(2)/(U(1) x U(l)) as for a moderate magnetic 
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field and spin degenerate levels. The difference in the prefactor a remains. One can extend this approach to other 
compact symmetric spaces with physical realizations, see Ref. ]3^,|33| for a complete classification. 

In addition to these gapless transversal modes there are massive longitudinal modes with Q 2 ^ 1, which for N 1, 
can be integrated out and the ASD thereby reduces to a functional integral over the transverse modes U. Now, 
the action of finite frequency to and spatial fluctuations of Q around the saddle point solution can be found by an 
expansion of the action F, Eq. (|60|). Inserting Q = UAU into Eq. ((6^), and performing the cyclic permutation of U 
under the trace Tr, yields, 

F = — - J dx < x | Trbx{Go 1 - U[H , U] + uoUMJ) | x >, (6) 

where 

G^=E-H + ^A. (7) 
Expansion to first order in the energy difference ui and to second order in the commutator U[Hq, U], yields, 

F[U] =-~u f dx<x\ TtGqeUMJ | x > 
+ l - J dx < x I TrG 0E U[H , U] | x > 

+ X - J dx < x | Tr(G 0E U[H , U]f | x > . (8) 

Note that [H , U] = -£(V 2 [7_) - f (Vt7)V - -^(r 3 AVC/ - C7t 3 AV). 

The first order term in U[H , U] vanishes for Gaussian white noise isotropic scattering. 

In general, in order to account for the ballistic motion of electrons in ballistic wires, or to account for different 
sources of randomness, a directional dependence of the matric U = Z7(x.n), where n = p/ | p |, has to be considered 
P^ , ^5| . However, for the geometries considered in this article, we have found that the form of the action derived below 
remains valid for diffusive as well as ballistic crossections, when the vector fields S as intorduced in Refs. p^j3q ], are 
integrated out. This will be presented in more detail in a separate article. 

Then, one can keep second order terms in V£7 and A, which turns out to be valid for the regime of weak disorder, 
/ ^> l/kp and for any magnetic field, Ib S> kp. Thus, one gets, using the saddle point equation, Eq. (||), 

r , n uj f dx . „ 
F M = -4Ajs-L TrAQ 

+ \J dx<x\ Tr(G 0E U(^-(WU)(W - ^Ar 3 ) + ^ [r 3 , C/AV])) 2 | x > . (9) 



Next, one can separate the physics on different length scales, noting that the physics of diffusion and localization 
is governed by spatial variations of U on length scales larger than the mean free path /. The smaller length scale 
physics, is then included in the correlation function of Green's functions, being related to the conductivity by the 
Kubo-Greenwood formula, 

fi q 2 

<Taf)(w) = FTT 2 " M(P{3 ~ 1 A p) < P I &0E I P' >< P' I G^E+u, I P >, (10) 

7T ij 1j TTh 

p.p' 

where p = §V. The remaining averaged correlators, involve products Gq E Gq E+uj and G^ E G^ E+LJ and are therefore by 
a factor H/(tE) smaller than the conductivity, and can be disregarded for small disorder Ti/t <C E. In the bulk of this 
article we are interested in the weak magnetic field limit, where uj c t <C 1, with the cyclotron frequency to c — qB/m. 
In this limit we can disregard the nondiagonal Hall conductivity and the explicit magnetic field dependence of the 
longitudinal conductivity. 

In order to insert the Kubo-Greenwood formula in the saddle point expansion of the nonlinear sigma model, it is 
convenient to rewrite the propagator in F as Gqe — \Gq E (1 + A) + \Gq E {\ — A). 

Then, we can use, that Tr[Y? a=1 E s =±( 1 + sK)U{V a U){l - sk)U{V a U)] = -Tr[(VQ) 2 ], and Tr[£ s =±(l + 
sA)U[t 3 , U]}(1 - sA)U[r 3 , U] = -Tr[[n, Q} 2 ]. 

Thereby we can rewrite Eq. (Q) as 
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7T LU f (fx , „ 

F M=-±AjsL TrAQ 

-\J dxTr[(VQ(x)) 2 < x | G «X(V - IgA)G^(V - JgA) | x > 

- -A 2 f dxTr[[r 3 ,Q(x)] 2 ] < x | G&AVG&AV | x > +c.c. (11) 
4 m J 



For wires of thickness W not exceeding the length scale Leu — Lc(fi = 2) = 2ttTivSDq 1 the variations of the field Q 
can be neglected in the transverse direction, and the action reduces to the one of a one- dimensional nonlinear sigma 
model. Using the Kubo formula, Eq. (10), this functional of Q thus simplifies, for lu c t -C 1, to, 



F = 77T2 (T ( UJ = °) W I dx(Tr(\7 x Q(x)) 2 - <A X *A X > ^Tr[T 3 ,Q(x)} 2 ). (12) 

io<r Jo ft 



The prefactor of the time reversal symmetry breaking term, the correlation function 

<A X »A X > = B 2 <yy> 

(13) 



(< x | G&AVG^AV | x > +c.c.) 



< x | G« £ (V - |?A)G^(V - f 9 A) | x > ' 

is increasing with the magnetic field B, suppressing modes with [Q, T3] ^ 0, the Cooperon modes, arising from the self 
interference of closed diffusion paths. Accordingly, the symmetry of the Q- fields is broken from Sp(2) / (Sp(l) x Sp(l)) 
to U{2)/{U{1) x [7(1)). 

In the next section it is shown that this prefactor is related to the magnetic phase shifting rate, and is evaluated 
for a disordered quantum wire. 



IV. THE MAGNETIC PHASE SHIFTING RATE 

It can be seen that the prefactor of the symmetry breaking term in Eq. ( |l2|) is proportional to the effective phase 
shifting rate 1/tb, governing the weak localization suppression by a magnetic field. To this end, one can use the 
supersymmetric version of the above nonlinear sigma model, obtained by substituting the matrix Q by supermatrices, 
and the trace over matrices Tr by the supertrace STr, but keeping all coefficients the same as in Eq. (|l2|). Then, the 
weak localization corrections to the conductivity can be calculated as outlined in Q], by an expansion of Q around 
the classical saddle point Q c = A. Thus, the magnetic phase shifting rate 1/tb can be identified as, 

Q 2 

1/tb = 4Z?^2 < A x • A x >, (14) 
n 

where the Einstein relation a = 2q 2 vD of the classical conductivity a to the classical diffusion constant D has been 
used. 



A. 2D wire with specular boundary conditions 



The general expression for the correlation function < y • y > , is found by inserting the momentum eigenstates of 
the wire and summing the correlation functions of Green's functions for Ib ^> W in Eq. (15). It is thus obtained to 
be given for a two dimensional wire of width W in momentum representation by, 



<W>= E kl{G* E (k x ,k y )G£ E {k x ,k' y )+c.c.) \<k y \y\k' y >\ 2 ±A x yG* E {k x ,k y )G£ E {k x ,k v ). 

(15) 

Here, G^ A (k x , k y ) = (E - h 2 (k 2 x + k 2 y )/(2m) ± i/(2r). 

Keeping all corrections for finite number of transverse channels N = k E W/ V and effective mean free path A = kpl, 
in the weak disorder limit E ^> U/r, we get for N 3> 1 the expression : 



G 



<yy> =W 2 (—K-—K 1 - 
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2tt 2 ' 



ir 2 N 2 



Kn 



N 



4 A 3 ^ s 



,ttN 



S— 1 



(16) 



where the definition of the constants i*Q is given in Appendix B . 

Its dependence on the mean free path parameter A = fcjrZ is shown in Fig. 1. 
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FIG. 1. The dependence of the correlation function < y • y > /W on the dimensionless mean free path A = k?l for 
N = 100 channels. For comparison, the line corresponding to a disorder independent phase shifting rate, approximately valid 
for N < A < iV 2 , is shown. 



Note that, although N 1 is required for the validity of the nonlinear sigma model, the equation ( |l6| ) is valid for 
arbitrary ratios of the width of the wire W and the mean free path I, since the motion remains diffusive along the 
wire axis on large length scales, even if I 3> W. 

For diffusive wire crossections, I < W, < y • y >— > y 2 = W 2 /12 which results exactly in the known result for the 

magnetic phase shifting rate 1/t b = 4D^y 2 B 2 [|],[|. 

The above derivation is more general, and applies for arbitary ratios of the wire thickness W and the mean free 
path /, as long as the magnetic length Ib is both larger than the width W and the elastic mean free path I, and for a 
large number of transverse channels N = fcf W/tt ^> 1. 

For ballistic wire crossections, I > W, Eq. (jig ) shows, that the effect of the magnetic field becomes weaker, as W/l 
decreases. This is a result of the flux cancellation effect, discussed in the limit of weak localization in Ref. The 
matrix element of the vector potential < k | A | k' > vanishes for k = k', since A = (—By, 0, 0) is antisymmetric in 
the coordinate perpendicular to the wire, y. Thus, elastic impurity scattering is needed to mix different momentum 
states and contribute finite matrix elements of the magnetic vector potential. 

One can check that Eq. (|l6j) is valid also in the weak disorder limit, by Taylor expanding the correlation function in 
l/(fcpZ), giving < yy >= ^-(A^ 3 /A 2 ), showing that it vanishes for A ^> A^ 2 , corresponding to h/r 
when the disorder does not mix transversal modes, like 1/A 2 , as seen in Fig. 1. 

In the intermediate regime, N < A, it had been argued in Ref. H|§, that 1/tb should be reduced by a factor linear 
in N/X resulting for a 2 dimensional wire with perpendicular magnetic field in a disorder independent expression 

1 1 W 3 v F 

— = — F (17) 

T B O l B 

where Ib = (h/ (qB)) 1 / 2 is the magnetic length. For specular boundary condition, as considered in this article, it was 
found numerically that C = 9.5 ||. Correspondingly, the function < y y > IW 2 should approach < y • y > /W 2 — > 
(tt/2C)N/\ or for N = 100, < y • y > /W 2 -> 16.5/A. The result Eq. ( [[J ) agrees indeed with this behaviour, 
in a regime N <C A <C A^ 2 , although the best fit gives a different prefactor 14.5, corresponding to C = 10.8. The 
analytical result shows, furthermore, that this behaviour is only an approximation and that there is a crossover to 
the perturbative regime, discussed above, where < y • y > /W 2 decays like ~ 1/A 2 , see Fig. 1. Note that this result 
is accurate upto corrections of order 1/N. 
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B. Parabolic Wire 



As long as the elastic scattering rate exceeds the cyclotron frequency, 1/r 3> u> c , or correspondingly I <C Icyc, 
where Icyc = kpl 2 B is the cyclotron path, determining the length scale on which ballistic paths start to bend due to 
the Lorentz force, the magnetic field dependence of the classical diffusion constant and the density of states can be 
neglected, being for a 2- dimensional wire D — and v(E) = m/ftirh 2 ), respectively. 

However, the cyclotron length can be small compared to the width of the wire, Icyc < W , while exceeding the 
elastic mean free path lc yc > I, when the crossection of the wire is diffusive, I < W . Thus, the localization length 
can depend sensitively on the ratio of these length scales, even in the weak magnetic field limit, where the density 
of states and classical conductivity are insensitive to the magnetic field. In order to study the crossover as function 
of the magnetic field, the dependence of the eigen functions on the magnetic field have to be taken into account, 
therefore. This regime is most conveniently studied for a parabolic wire, having a harmonic confinement, 



having the energy eigen values 



H ° = 2m" (P ~ qA)2 + l muJ ° y2 > ( 18 ) 



E n , k = huj cB (n + 1/2) + ^h 2 k 2 , (19) 

2m* 



where the effective mass is m* = mui 2 s /uj 2 , and the effective frequency is ui e g — (cj^+u; 2 ,) 1 ' 2 , where lob = qB/m is the 
cyclotron frequency. The spatial center of the electron eigenstates are shifted by the guding center yf. = khus/ ( m ^ 2 ff ) ■ 
Thus, the width of the wire is at constant Fermi energy Ep dependent on the magnetic field B. Defining the width 
of the wire W at fixed Fermi energy as W 2 = max(< n, k | y 2 | n, k >) with E n ,k = E-p, one finds for the parabolic 
wire: 

E V >< 2 



W 2 (B) = Z c ynax(2-^^ + („ + 1/2)(1 - %). (20) 



For large magnetic field, lob S> <-^o, this approaches exactly twice the value at zero magnetic field, and thus, 

W(lo c > uJo) = V2W{0) = (2E F /(huj )) 1 / 2 l . (21) 

Thus, the wire width is a slowly vaying function of the paramter lu c /luq = W(B = 0)/7 cyc . 

The presence of impurities smoothens this function further, and we can thus assume the width to be practically 
magnetic field independent: 



W = y/2E F /m e /u; . (22) 

This allows us to study the various regimes of interest as a function of the wire width W, the magnetic length Zb and 
the average mean free path I = (2E/m) 1 / 2 T. 

Naturally, the classical conductivity in such a wire is anisotropic. We find that 



1 



J c« T 



q 2 Tn e /m, (23) 



and 



'vv 



T q 2 Tn e /m, (24) 



where n e — (2/3n)(m e E/?i 2 u;o) is the average electron density in the wire, which is taken to be approximately 
independent of the magnetic field. Since we consider magnetic fields where wqt <C 1, the classical conductivity is 
magnetic field independent, a X x — q 2 Tn e /m, and cr yy — cr X x/(l +lo 2 t 2 ). 

Thus, the condition that the localization is governed by the one-dimensional nonlinear sigma model is changed to 
Leu /(I + w o t2 ) > W. With loqt = l/W follows that the one dimensional localization condition requires, / < 2NW, 
in the weak disorder regime, kpl 1. 

Rederiving the nonlinear sigm model in the representation of a clean parabolic wire, using the definition of the 
correlation fucntion, Eq. dig), where teh sum over transverse momenta is substituted by the sum over the band index, 
n, k y —>■ n, we find the result, 
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2 1 ui 2 2 1 W 2 

<y ' y> - w2 ^TT^ + 3 4 } = W2 5 ( TT7W + V' (25) 

Note that, since uj^t 2 — l 2 /(W 2 ), the ballistic crossection limit I > W, coincides for the parabolic wire with the 
clean wire limit, where transversal modes are not mixed by the disorder h/r < Hluq. Thus the flux cancellation effect 
leads in the parabolic wire to a supppression of the phase shifting rate by a factor W 2 /l 2 as found for the wire with 
specular boundaries in the clean wire limit as seen in the previous subsection. 

Thus, it is not surprising that the behaviour of the magnetic phase shifting rate, as known from weak localization 
corrections for a wire with ballistic crossection, W > I, and hard wall boundary conditions, is not reproduced when 
considering a parabolic wire. In the former case, there is a regime, W 2 < l 2 B < Wl, implying Is < I, where the 
magnetic phase shifting rate is given by 

1 W 2 W 3 v F 



where C*2 = 24/5. This is smaller than expected from Eq. (|17|), and is not obtained for the parabolic wire. 

Instead, we find that there is a regime, where the magnetic field sensitivity of localization becomes stronger, when 
the cyclotron length i cyc , becomes comparable to the width of the wire W . When I < l cyc < W the magnetic phase 
shifting rate is found to increase with the magnetic field like B A , 



1 24 q 2 9 W 4 
— = —D-, T '~ 
tb 5 h 



= ~ D h BZ — - (27) 

eye 



When the magnetic field becomes so strong that the cyclotron length Z cyc , becomes comparable or smaller than the 
mean free path I, or lo c t > I, the diffusion constant and the density of states become functions of the magnetic field. 
Then, the spatial modes of the nonlinear sigma model perpendicular to the wire can become soft and contribute to 
the functional integral, and thus, the nonlinear sigma model becomes effectively two dimensional. 

In this limit, a quantum Hall wire, the approach used in this article can yield qualitative information on the location 
and size of localized states in a quantum Hall system [E9|, and will be reconsidered in a forthcoming work. 



V. MAGNETOLOCALIZATION IN DISORDERED QUANTUM WIRES 

It is known that the localization length depends on the global symmetry of the wire ||: L c = flnhvSDo, where 
(3 = 1,2,4, corresponding to no magnetic field, finite magnetic field, and strong spin- orbit scattering or magnetic 
impurities, respectively. v(E) is the electronic density of states in the wire |l],|2|]. Dq is the classical diffusion constant 
of the electrons in the wire, and S its crossection. This result was obtained by calculating the spatial decay of the 
density correlation function for wires whose thickness exceeds the mean free path I. 

Here, we use an extension of a recent nonperturbative calculation, to obtain the localization length as a function 
of the magnetic field, using the fact that the ASD shows a crossover from an oscillating behaviour, decaying with 
a power law [p|JlCj], typical for Wigner- Dyson energy level statistics to a gaussian decaying function, when the 
length of the wire is increased beyond the localization length |Tq| , as seen in other measures of correlations in the 
discrete energy level spectrum of a phase coherent disordered electron system P, ^6| , pTp^JT4] ] . 

Taking the representation of the ASD derived above, Eq. (f59|), 

C(u>) = J JJdQWexp(-F[Q]), (28) 

where the action Eq. ( |l2| ) can be rewritten conveniently in terms of the diffusion length, an electron would diffuse 
classically in the magnetic phase shifting time tb, Lb = \/Dtb'- 

r-L 



F[Q]=a^L CU J o dxTr 



{V x Q{x)) 2 - -^t[QM 2 



+ J ^TrA 3 Q(x). (29) 

where Leu = Lc{f3 = 2) = 2irUvSD is the localization length in the wire in a moderately strong magnetic field ||. 

In the limit when Lb < Lc, a moderately strong magnetic field, Q is reduced to a 2 x 2- matrix by the broken time 
reversal symmetry. This reduces the space ofQ to U(2)/(U(1) x U(l)). 







For ix> I A < Lcu/L, corresponding to ui < Eq, where Eq — 2-kD/L 2 is the Thouless energy scale of classically free 
diffusion through the wire of length L, the spatial variation of Q can be neglected and one retains the same ASD as 
for random matrices of orthogonal or unitary symmetry, respectively ||,[l0| . 

Increasing the length of the wire L, a crossover in the autocorrelation function can be seen as the wire exceeds the 
length scale L c . 

In order to study quantum localization along the wire, the function C(cu) should be thus considered as a function 
of the finite length L of the wire and spatial variations of Q along the wire have to be considered, as described by the 
one dimensional nonlinear sigma model derived above. 

The impurity averaged ASD can to this end be written as a partition function [ p9[ 

C{uj) = Trexp(-LiJ [Q]), (30) 

where H is an effective Hamiltonian of matrices Q on a compact manifold, determined by the symmetries of the 
Hamiltonian H of disordered electrons. Thus, the problem reduces to the one of finding the spectrum of the effective 
Hamiltonian H . 

We can derive the corresponding Hamiltonian H by means of the transfer matrix method, reducing the one- 
dimensional integral over matrix field Q, Eq. fl59|), to a single functional integral. Thus, the ASD is obtained in the 
simple form of Eq. ( |30| ) , with the effective Hamiltonian 

H{u = 0) = (-4A* - ±-X 2 Tr Q [Q, r 3 ] 2 ). (31) 

aL C u 16 

Aq is that part of the Laplacian on the symmetric space, which does not commute with Tr[A3<5]. The time reversal 
symmetry breaking due to the external magnetic field is governed by the parameter X = aLcu I '(2£s). 

The problem is now equivalent to a particle with "mass" (a/8)Lcu(E) moving on the symmetric space of Q in a 
harmonic potential with "frequency" 1/(2Lb), and in an external field za(7r/4)w/(LA), in "time" x, the coordinate 
along the wire. To find the ASD as a function of lu and the length of the wire L, one can do a Fourier analysis in 
terms of the spectrum and eigenfunctions of the effective Hamiltonian at zero frequency, H(ui = 0) Jjfl. 

There is a finite gap Eq between the ground state energy and the energy of the next excited state of H{lo = 0). 
For a long wire, LEq ^> 1, the ASD becomes, C(lo) = exp(— const. Loj 2 / Eq) , where both const.u 2 =|< 
H(u)) — H(0) | 1 >| 2 , and the gap between the ground state and the first excited state, E G = E x — E Q do depend on 
the symmetry of the Hamiltonian H . This exponential decay with Loj 2 is typical for a a spectrum of localized states 
p9| . In the other limit LEq <C 1, all modes of H do contribute to the trace in the partition function Eq. (|3(]) with 
equal weight, yielding the correlation function of a spectrum of extended states Jig] . Thus, the crossover length is 
entirely determined by the gap Eg, through £ c = 1/Eq, and can be identified with an averaged localization length. 

In order to derive the eigenvalues of the effective Hamiltonian at zero frequency, H(u> = 0), we need to introduce a 
representation of the matrix Q and evaluate the Laplacian in its parameters. This is done in Appendix C. 

Without magnetic field, B — 0, the Laplacian is obtained to be 

A£ = Xc (1 - X 2 c )d Xc + 2^ T ^5 Ac 



A 

^LdxJl-^Ddxv, (32) 
A c 

where Ac,d € [—1, 1]. Its ground state is 1 and the first excited state is Xc^d- Thus, the gap is 

E G (B = 0) = 16/Lcu- (33) 

For moderate magnetic field, with the condition Lcu{< U'V >) 1 ^ 2 B ^> (f> — h/q, all degrees of freedom arising from 
time reversal invariance are frozen out, due to the term TtqIQ^t^] 2 = 16(Ap — 1) which fixes A^ = 1. Then, the 
Laplacian reduces to 

A%=d XD {l-\ 2 D )d XD . (34) 

Its eigenfunctions are the Legendre polynomials. There is a gap above the isotropic ground state of magnitude 

E G (X > 1) = 8/ L cu - (35) 

For moderate magnetic impurity scattering, exceeding the local level spacing, 1/t s > Ac, a = 2, and the Laplacian 
is given by Eq.(p3). 
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Thus, due to a = 2, the gap is reduced to Eq(1/ts > Ac) = 4/ Leu- For moderately strong spin- orbit scattering 
l/ T SO > Ac, the Laplace operator is 



A* = J2 0a,(1-A?)0* 

Z=l,2 



(36) 



where Ai^ G [—1, 1]. The ground state is -00 = 1, the first excited state is doubly degenerate, ipn = Ai, ipi2 = A 2 . 
Thus, the gap is the same as for magnetic impurities, 



E g (1/t S o > A c ) =4/L C u- 
An external magnetic field lifts this degeneracy but does not change the gap. 

TABLE I. Relation between symmetry of the Hamiltonian and the gap of the quasi-lD- NLSM 



(37) 



Class 


I Symmetry 


Symmetric Space 


Cartan class 


Gap Eq 


Ordinary 


T R 


S R 


Sp(2)/( Sp(l) x Sp(l)) 


CII 


W/Lcu 


Ordinary 


No T R 


S R 


t/(2)/(f/(l) x [7(1)) ( Sphere ) 


AIII 


8/ Leu 


Ordinary 


T R 


No S R 


0(4)/(0(2) x 0(2)) 


BDI 


4:/L C u 


Ordinary 


No T R 


No S R 


U(2)/(U(1) x [/(!)) 


AIII 


4/Lcu 



Thus, using the crossover in energy level statistics as the definition of a localization length as above, we get in a 
quasi- 1 -dim. wire, 



& = l/E B {p) = (l/16)pLcu, 



(38) 



where (3 = 1,2,4 corresponding to no magnetic field, finite magnetic field, and strong spin- orbit scattering or 
magnetic impurities, respectively. Comparing with the known equation for the localization length, L Cl we find that 
the dependence of the ratios (3 on the symmetry are in perfect agreement with the result as obtained from the spatial 
decay of the density- density- correlation function ||, while it defers by the overall constant 1/8. 








L/2 



a) 



b) 



x 

FIG. 2. Schematic visualization of the energy level spectrum of localized states in a) a disordered quantum wire of length 
L, when divided into two parts, b) for the same wire when both parts are connected and the eigenstates are hybridized. 
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This relation can be proven directly. The ASD at zero frequency C(0)l of the wire of length L, becomes, when the 
wire is divided into two parts, C(Q)^i 2 . For L — > oo, we find that the relative difference is: 

f(L) = -~^j£-l = 2pe X p(-LE G /2), (39) 

exponentially decaying with the length L. Here p is the degeneracy of the first excited state of H(u — 0). f(L) 
can be estimated, following an argument by Mott p7fl : When the two halves of the wire get connected, see Fig. 
(2) , the Eigenstates of the two separate halves become hybridized and the Eigenenergy of a state ip n is changed by 
±Ac exp(— 2x n /Lc). x n is random, depending on the position of an eigenstate with closest energy in the other half 
of the wire. Thus, averaging over x n gives: 

/(L)~+exp(-4L/L c ). (40) 

Comparison with Eq. ( |39| ) yields indeed 1 / Lc = SEq ■ 

It is thus a remarkable fact that this length scale, defined as the crossover length of the spectral autocorrelation 
function, and related to the excitation gap of the compact nonlinear sigma model, has exactly the same symmetry 
dependence as the localization length, defined through the exponential decay of the spatial density correlation function, 
found in Ref. || . This is especially surprising, since the nonperturbative derivation of the disorder average of the 
quantity, < p(r, t)p(r', t') > — < p(r) 2 >, necessitates the use of the supersymmetry method, resulting in a nonlinear 
sigma model of supermatrices, having in addition to a compact sector, the one considered here, a non compact sector, 
where the matrix is parametrized on a semi infinite interval. The full supersymmetry allows furthermore rotations 
between this compact and noncompact sector which are parametrized by Grassmann numbers £, having the property 
£ 2 = 0. Apart from this increase of the manifold of the matrix fields Q to the supersymmetric space, the structure 
of the theory is equivalent. Especially, the free energy of the supersymmetric nonlinear sigma model, has exactly the 
same form as Eq. (|29|), replacing Q by supermatrices and the Trace over Q, by a supertrace STr, giving the opposite 
sign to the noncompact sector |2| . 

Studying localization in a wire with this supersymmetric nonlinear sigma model, the transfer matrix method yields 
an effective Hamiltonian of supermatrices Q, of the same form as Eq. (^), where the Laplacian is now defined on 
the respective supersymmetric manifold. In full analogy, the spectrum of H determines accordingly the properties 
of a disordered quantum wire, and has been derived in Ref. 3(| for the pure ensembles. The partition function 
Z = STrexp(-LH) is a generating function of spectral correlation functions 38,17]], In order to derive spatial 
correlation functions like the density correlation function, in addition, the Eigen functions of the respective diffusion 
equation on the supersymmetric manifold, 

(-d x +H(Q))ip(x;Q) = 0, (41) 

have to be found || . In that way, a formula for the conductance of a finite disordered wire attached to two leads at a 
distance L, has been derived |56|, see also Ref 0. In the limit of a wire which is perfectly coupled to the leads, that 
formula for the average conductance simplifies to 

< 9 >= [ dfi(k)E(h) exp(~E(!i). (42) 
2a J 16 

Where E{li) are the eigenvalues of the supersymmetric Hamitlonian H (u> — 0) and dp(h) the corresponding integration 
measure, of the discrete and continous eigenvalues of the angular momentum operator on the compact and noncompact 
sector, respectively. They were found to be given for B = by 

E(k) = 0, 4/L cu 2(e 2 + l)A/L CU {l 2 + 4+4 + 1), (43) 

where I = 3,5, and e > 0, t\ > 0, 62 > 0. 

For time reversal symmetry broken wires X > 1 the eigenvalues were found to be, 

E(k)=0,— —M 2 + e 2 ), (44) 
otLcu) 

where Z = 1,3, 5, and e > 0. 

If spin symmetry is broken, but time reversal symmetry conserved, in the presence of spin orbit scattering, the 
eigenvalues were found to be, 
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E(k) = 0, ^-(2(/ - I) 2 , + l\ + ? 2 - 1), (45) 

where I = 3,5, U = 1, 3, 5, = 1,2 and e > 0. 

In that case it can be seen that for a distance between the leads much exceeding the localizaion length, L 3> Leu, 
the conductance decays exponentiallly, and that this is entirely determined by the compact gap Eq, between the 
lowest angular momentum eigenstates of the compact sector. The integration over the continous eigenvalues of the 
noncompact sector, leads only to a prefactor, decaying as a power of the length, ~ 1/L 3 / 2 . Indeed, the gap between 
the ground state value E = and the first excited state is seen from Eqs. (H@||), to be E G = 8/L cu for B = 0, 
Eq = 4/ Leu for X > 1, E G = 2/ Leu f° r magnetic impurity scattering, a = 2, and Eg = 2/ Leu for moderate 
spin- orbit scattering, coinciding with the symmetry dependence of the compact gap derived above. However, that 
coincidence might appear as mere chance, since in fact, the Laplacian of the supersymmetric matrix Q can not be 
written as a sum of the one of the respective compact nonlinear sigma model, Eqs. ( ^2| , p4|j3^ ), because the metric 
tensor g on the supersymmetric space contains mixed factors of compact and noncompact parameters. Therefore, the 
discrete eigenvalues of — Aq, are not the Eigenvalues of the square of the angular momentum on a compact sphere 
p6| . Only, in the limit of infinite noncompact parameters does one recover the respective Laplacian on the compact 
symmetric space, Eqs. ( |3^ , |34] , |36| ). 

Thus, having shown that the ASD yields the correct symmetry dependence of the localization length, we can now 
use this approach to get an analytical solution for the crossover behaviour of the localization length and the local 
level spacing as a magnetic field is turned on, and there is no spin- orbit scattering. While a self consistent approach 
p9[ , a semiclassical analysis |4(J] and numerical studies |^l|,Q showed a continous increase of the localization length, 
an analytical result indicated that both limiting localization lengths L c ((3 — 1) and L c (/3 = 2) are present in the 
crossover regime and that there is no single parameter scaling. This is explained by arguing that the far tails of the 
wavefunctions do cover a large enough area to have fully broken time reversal symmetry, decaying with the length 
scale L c (f3 = 2) even if the magnetic field is too weak to affect the properties of the bulk of the wavefunction, which 
does decay at smaller length scales with the shorter localization length L c (f3 = 1), corresponding to the time reversal 
symmetric case. The quantity studied there is the imprurity averaged correlation function of local wavefunction 
amplitudes and its momenta at a fixed energy e: Y(e) =< J2 a I ^o(O) | 2 | V'aM | 2 ^( e — e a) >■ It is averaged 
over a distribution of eigenfunctions in different impurity representations. Thus, each eigenfunction could decay 
exponentially with a single localization length, but having a distribution which has two maxima, at L c {j3 — 1) and 
L c (f3 — 2), whose weight is a function of the magnetic field in the crossover regime. While the distribution function 
of ln(| i/j a (0) | 2 | ip a {r) | 2 ) is known to be Gaussian in both limiting cases of conserved, and fully broken time reversal 
symmetry, centered around the value r/Lc(/3), (3 — 1,2, respectively, it is not yet known in this crossover regime, 
however The average value of moments, | ip a (Q) \ k \ ipa{r) \ k , is decaying more slowly than its typical value, and 
does not depend on the order of the moment, k. This was taken as a proof that moments are determined b y st ates with 
anomalously large localization lengths of the order of the system size |ll|. Therefore, the result of Ref. can be a 
property of such rare states with anomalously large localization length, and it remains to see, if the full distribution 
function scales with two lengths L c ((3), f3 = 1, 2, or a single one, changing continously with the magnetic field, L C (B). 

While we cannot resolve this question by calculating a spectral autocorrelation function like the ASD, this is another 
motivation to see if the energy level statistics is governed by a single parameter as the magnetic field is varied. 

The effective Hamiltonian for moderate magnetic fields is found, without spin dependent scattering, a = 1, using 
Tr[Q, r 3 ] 2 = 16(1 - X 2 C ) to be given by: 

^-7^(-4A0 + A 2 (l-A 2 7 )), (46) 
Leu 

where the Laplacian is Eq.(j3^) and X = Leu I '(2Lb)- 

In the limit X — > the ground state and first excited state approach 1, Ac Ad, respectively. In the limit X ^> 1, X c 
becomes fixed to 1. Thus, the Ansatz V'o(^c) ~ exp(A X 2 (l — X c )), and tpi(Xe, Xo) ~ A C A D exp( J 4 1 A 2 (l — A^)), 
where Ao < 0, A\ < are negative constants, solves Hip — Eip to first order in z = X 2 (l — X c ). One finds that the two 
lowest magnetic field dependent eigenvalues are Eq = 4/Lcu(—5 + v25 + X 2 ), and E\ = 4/Lcu{—3 + V49 + X 2 ), 
and the Eigenfunctions are given as above with Aq = —LcuEq/(16X 2 ), and A\ = (1 — LeuEi/16)/X 2 , yielding the 
right limits for X — > and X ^> 1, respectively. Thus, there is a magnetic field dependent gap Eg — E\ — Eq of 
magnitude: 



E G (X) = 4(2 + ^49 + X 2 - ^25 + X 2 )/L CU - (47) 
This solution is valid in both the limits X <§C 1 and I>1, interpolating the region Iral, 
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With the magnetic diffusion length Lb = (Dtb) 1 ^ 2 , and the magnetic phase shifting rate, as given by Eq. (|l4|), we 
obtain: 

X = L CU /(2L B ) = L cu j:V< V • y >B, (48) 

which is y/< y • y >/W times the number of flux quanta penetrating a localization area LcuW. 

From Eq. ( f47| ) follows that the magnetic change of the localization length is SLc{B) ~ B 2 for small and ~ l/B at 
large magnetic fields, which agrees with the result of the selfconsistent method as obtained by Bouchaud p9]. 



VI. RESISTANCE OF DISORDERED QUANTUM WIRES 

In the limit of zero temperature, T = 0, the resistivity of a disordered quantum wire, having only localized 
states at the Fermi energy, is infinite. For finite temperature, T > 0, in the strong localization regime k-gT < Ac, the 
mechanism of conduction is hopping of electrons between localized states. Then, the resistivity increases exponentially 
with temperature. According to the resistor network model each pair of localized states i and j is linked by 

a resistance Rij : 

% =exp(^ + -^-) (49) 
L c k B T 

where = jr* — rj\ and ejj = (|ej — /Lt| + |e 3 — /i| + |ej — ej\)/2ksT (r^ and ej are the position and energy of the state 
i, fi being the Fermi energy). Because of the exponential dependence of R on and ey, percolation theory methods 
can be applied ]4^j4^|. In 2-D and 3-D systems, the dependence of R on temperature T shows a crossover from an 
activated behaviour to the variable range hopping (VRH) regime. In this regime the temperature is so low that the 
typical resistances between neighbouring states are large because of the second term in Eq. ( f49|) . Therefore electrons 
tunnel to distant states whose energies are close to the Fermi level. If we neglect electron-electron interactions the 
resistivity is described by Mott's law |49| , ^| : 

R(T) = R o exp[( 7 T /T) 1 /^] (50) 

where d is the dimensionality of the system, 7 a numerical coefficient which depends on d, To = and Vd is the 

dimens dependent density of states. However, in the quasi- 1-D case and for sufficiently long wires the variable range 
hopping result, Eq. (pp[), cannot used due to the presence of exponentially rare segments inside which all the localized 
states have energies far from the Fermi level |5(] -|5^]. These large resistance segments (LRS) do not strongly affect the 
resistivity of 2-D and 3-D systems because they can be circumvented by the current lines. In 1-D this is not possible 
and the total resistance of a wire is given by the sum of the resistances of all the LRS's. This sum yields an activated 
type dependence of R on T ]5j| for infinite wires: 

i? = i?0 ^ ( Y ) ^ exp(To/2T) ' (51) 



where fcsJo = l/vL c = A c coincides with the local level spacing, and L is the length of the wire. Eq. (51) is valid 
provided that the number of optimal LRS's (i.e. those LRS's which give the largest contribution to R mW) within the 
length of the sample) is large. Bur for a finite wire length this condition fails to be fulfilled at very low t emp erature 
T, and the resistance of the chain is determined by smaller LRS's; in this regime Eq. © is replaced by HJH): 



/^« xi»[\/2^1og(-^(^)hogi(£))] I (52) 



which is valid below a temperature 



T 

Tl = 2 h\(L/L c ) ' (53) 



approaching Mott's law, Eq. ( p30| ) at lower temperatures T < T±. 

So far, electron-electron interactions have not been taken into account. This approximation is valid if the Coulomb 
interaction is screened over distances of the order of the hopping length, as by a metal gate electrode deposited 
on top of the wires at a distance smaller than the typical hopping lengths. When this is not the case, long range 



electron-electron interactions affect both the density of states and the resistance of the samples 53 5 
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VII. COMPARISON WITH EXPERIMENTAL RESULTS 



The magnetic field dependent activation energy was measured recently in transport experiments of Si 5- doped Ga 
As quantum wires [0. As an example, we discuss here the sample 5 of Ref. 0, with a width W = .2fim, a localization 
length Leo = -61/iTO a length L = 40/j,m, and N = 30 channels. 

The activation energy coincides with the local level spacing fee To = A c = l/(vWL c ) and is estimated for sample 5 
to be T = MK. 

Thus, according to the theory outlined in the previous section, there is an activated reistance in an order of 
magnitude temperature range T\ = .04A' < T < Tq = .3AK , allowing in good approximation the direct measurement 
of the magnetic field dependent activation energy A C (B), and thus the magnetic field dependence of the localization 
length L C (B). 

The ratio of the cyclotron frequency and the elastic scattering rate, loqt — l/(kFl%) -C 1, is small in the whole 
range of magnetic fields considered there, so that the classical conductance would be magnetic field independent, 
(7 = ne 2 r/m(l + lo c t 2 )~ 1 m ne 2 T jra. 

The mean free path I ~ .02^ is small compared to the width of the sample W = .2fim. The magnetic length is 
Ib = ■026fim(B/T)^ 1 / 2 . Thus, while luc t <C F the magnetic length becomes smaller than the width of the sample 
at magnetic fields B > .0165T. 

The experimental magnetic field dependence of the ratio of activation energies is shown in Fig. (3) together with 
the theoretical curve for the ratio of local energy level spacings Ac(B)/Ac(0) — Eg(B)/Eg(0), as derived above, 
Eq. ( 36 ), using for the magnetic phase shifting rate the results for a 2-dimensional wire with specular boundary 
conditions, Eq. (14), and, for comparison, the one derived for a parabolic wire, Eq. (|25|). 

There is a quantitative discrepancy between the best fit X = .036B/G, and X = 2n4>/<po, 4> = ^oHLcviu 2 )^ 1 ^' , 
when using the analytical formula Eq. ( 14). With the experimental parameters a — I, Leo — -61/im, width 
W = .2fim of sample 5 in Ref. Q and y 2 = W 2 /12 for a 2- dimensional wire, it yields rather X = .OWB/G. We note 
that smooth confinement can give y 2 > W 2 /12. A similar discrepancy was observed between W as obtained from the 
sample resistance and estimated from the analysis of the weak localization magnetoresistance, which also depends on 

?!• 

We note that the agreement, when using the experimental parameters, for the parabolic wire, is better. The 
cyclotron length l cyc = kpl^ = .32/ (B/T)fim, is found to be larger than the mean free path I for B < 15T and larger 
than the wire width for B < 1.5T. We find for the parabolic wire: X = -024^99 + 1.33 10~ 8 (B /G) 2 ) 1/2 B /G. The 
enhancement of the magnetic phase shifting rate in a parabolic wire, Eq. (Eq), is thus too weak to be seen at the 
magnetic fields used in the experiment, B < .2T, as shown in Fig. (3), and seems thus not to be the origin of the 
increase in the decay of the activation gap, at about .IT . 

An extension of the derivation given in section IV to include a dependence of the eigenfunctions on the magnetic 
field also for a 2-dimensional wire with specular boundary conditions has to be done, in order to make the comparison 
with the experiment more quantitative, and conclude from the magnctolocalization on the form of the confinement 
potential in these Si- S- doped Ga As quantum wires. But, our results may indicate that the harmonic confinement 
model of the parabolic wire is a better description of the wires in sample 5. 



c 
g 
S 



Parabolic Wire 
sample 5, T=.3 K 
X=.01 
X = .036 




10000 



FIG. 3. The activation gap ratio To(H) /To(0) as a function of the magnetic field B in G of sample 5 measured at temperature 
T = .3K as reported in Ref. 1, together with the theoretical curves for a parabolic wire, using the parameters of sample 5. and 
a 2D wire with specular boundary conditions for a best fit value X = .036B/G, and the value obtained from the experimental 
parameters, X — .010B/G. 
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VIII. SUMMARY AND OPEN PROBLEMS 



A formula for the magnetic phase shifting rate has been derived, which allows its calulcation for arbitrary wire 
geometries and ratios of the elastic mean free path, the wire width, and the magnetic length. 

For a quantum wire with specular boundary conditions and harmonic confinement this formula has been evaluated 
explicitly, and compared with previous analytical and numerical results for the magnetic phase shifting rate. 

The localization length is derived as the crossover length scale from correlated to uncorrelated energy level statistics, 
as studied with the autocorrelation function of spectral determinants. It is shown that its symmetry dependence 
coincides exactly with the localization length as defined by the exponential decay of the averaged two-terminal 
conductance and derived with the supersymmetry method. 

Therefore, the ASD can be used to get analytical information on the magnetic field dependence of the localization 
length, which is shown to be governed by the magnetic phase shifting rate, and thus strongly dependent on the 
geometry of the wire and the ratios of the elastic mean free path, the wire width, and the magnetic length. 

A comparison with the magnetic field dependence of the activation gap, as observed in low temperature resistance 
measurements in Si 8- doped Ga As wires, indicates, that the electrons move in a potential which is closer to a 
harmonic than a hard wall confinement. 

Enhancement of the sensitivity of the localization to a magnetic field is found analytically, when the cyclotron 
length is comparable with its width. The physical reason for this enhancement is found to be the magnetic field 
dependent shift of the guiding centers of the electronic eigenstates in the quantum wire, even at moderate magnetic 
fields, when the classical conductivity is still independent of the magnetic field. 

It remains to extend the derivation to include random surface scattering EJ and the effect of correlated, smooth 
disorder ]57{ , in order to allow for a more quantitative comparison with the experiment. Both effects necessitate a new 
derivation of the nonlinear sigma model, which allows for a directional dependence of the matrix field Q. This has 
been recently introduced for a system with broken time reversal symmetry in the study of localization in correlated 
disorder j34|, and the spectral statistics of quantum billards with surface scattering In both cases one is lead to 
a nonlinear sigma model, where variations of the matrix Q on ballistic length scales are taken into account [p8f-|60[| . 
The application of this approach to the magnetolocalization in disordered quantum wires will be presented in a future 
publication. 



The authors gratefully acknowledge, usefull discussions with Isa Zarekeshev, and Konstantin Efetov, thank Yuri 
Khavin for providing his data and usefull communications, and Bernhard Kramer for stimulating support and critical 
reading of the manuscript. 



Here, the derivation for spinless case, a = 1, is given in detail. We use for compactness the vectors of anticommuting 
variables, 



Note that ip = (Cip) T , where the matrix C interchanges the Grassmann fields with their conjugate one, and has thus 
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Thus, the ASD is written as 



1G 



c(w)= J n^(x) 



exp(-i J dx$(x)(E + ^wA -H - 7(x))V(x)), (55) 

Here, the diagonal Pauli matrix A = ^ ^ ^ J has been introduced for compactness, its diagonal elements pro- 
jecting on the respective spectral determinant of the ASD. The kinetic Hamiltonian becomes a matrix 

H Q = (p-qT 3 A) 2 /2/m, (56) 

where the diagonal Pauli matrix T3 had to be introduced since each vector has elements of the Grassmann field and 
the time reversed one, and the diamagnetic term pA/m in the Hamiltonian changes sign, as p — * — p, breaking the 
time reversal invariance. 

To summarize the notation, here, and in the following, Aj are the Pauli matrices in the subbasis of the left and 
the right spectral determinant, the ones in the subbasis spanned by time reversal and <ii the ones in the subspace 
spanned by the spinor, for i = 1,2,3. 

Note that a global transformation of the Grassmann vectors rp =4- ?/> = Atp does leave the functional integral for 
uj = invariant, as long as A + A = 1, and A +T C = CA, restricting the matrices A to be symplectic ones, being 
elements of Sp(2), commuting with the antisymmetric matrix C . A finite frequency breaks this symmetry group, and 
only symplectic transformations of each field of a single spectral determinant separately, Sp(l) x Sp(l), do leave the 
functional integral invariant. 

Now, the averaging over the disorder potential can be done, integrating Eq. (|5^) over the Gaussian distribution 
function of the random potential V. 

Thus, the averaged ASD is found to be given by a functional integral over interacting Grassman fields ip, 



exp(-i J d^{x)(E + ^cjA-p 2 /2/m)^)) 

exp(--i- — SL f dxTr(^(x) x ^(x)) 2 ). (57) 

l07T T J 

Now, the resulting ^-interaction term can be decoupled by introducing another Gaussian integral over an auxilliary 
field. Clearly, the field should not be a scalar, otherwise we would simply reintroduce the Gaussian integral over the 
random potential V. Rather, in order to go a step forward, the auxilliary field should capture the full symmetry of 
the autocorrelation function. Therefore, the Gaussian integral is chosen to be over a 4 by 4 matrix Qixi, which is 
itself an element of the respective symmetric space, as the matrix A which leaves the functional integral invariant. 
Thus, allowing for a spatial dependence of Q, one can decouple the interaction term: 

\2\ 



exp(--L — SL f dx(V(x) x ^(x)) 2 

1D7T T J 

J nd<24x4(x)exp(-7r^ J ^-TrQ 4xi (x) 



+ i^J dxTrQ 4x4 (x)^(x) x V^(x)). (58) 

Anticipating, however, that the functional integral over the matrices Q cannot be performed exactly, but rather only 
an integral over slowly varying modes around a saddle point solution, it is necessary to separate fast and slowly varying 
modes already before the decoupling of the interaction term Eq. (^8|) |22| . It turns out that there are two equivalent 
slowly varying interaction terms, corresponding to diffusion, and one arrives finally after a Gaussian decoupling to a, 
by a factor 1/2, shallower nonlinear coupling TrQ 2 0. 

Next, one can perform the Gaussian integral over the Grassmann vectors V'(x) and one obtains for the ASD, 
rescaling Q4X4 — » 2t/HQ4 X 4, the representation: 

= /n d< 34x4(x)exp(-F[Q]), (59) 



17 



with 



+ dx < x | Tr\n(G(x,p) | x >, (60) 



where 



G(x,p) = 1/(~ W A - (P 2 ^ A)2 + ? ^Q 4 x 4 (i)). (61) 

is the propagator matrix. We used the operator notation x, in order to stress that the terms in the inverse propagator, 
do not commute with each other. 

APPENDIX B 

For a clean wire with hard wall boundaries, the transversal Eigcn modes arc for — W/2 < y < W/2, < k y | y >= 
cosk y y for k y = irs/W, s being an odd integer, and < k y \ y >= smk y y for k y = tts/W, s being an even integer, one 
obtains: 

^ v ^ i= w^v^-v^ )2 ' (62) 

when k y — ns/W , and k' y = us' /W , s being even, and s' odd, or vice versa. Then, the sum over k' y in Eq. ( [l5] ) can 
be performed by use of the Matsubara trick, for s even, and odd integers, separately. The remaining sum over k x , k y 
can be transformed as l/(W£)Efc x k = J dev(e) J noting that the unit vector e/. can point only in discrete 
directions. Thus, while in 2 dimensions J ^ = Jq" = 4/ (2tt) Jq dyl/{\ — y 2 )^l/2) 1 for finite number of transverse 
channels N = k F W/n there is a sum, J = 2/(nN) J2 S>0 ~ s 2 /N 2 )^l/2). Thus, k y = ns/W = kps/N 
and k x — kj?(l — s 2 /N 2 )(\/2). Performing finally for E ^> %/t the integral over e, one arrives with some patience 

n A ~ r_o//„ A n^ . f, ~ _ i//,AAr» . /1 s 2 /„2 



at Eq. f§, where A T = 2/ jnN) £^ =0 V 1 " #n * = VMO £"=x V 1 " W> K ^ = 2 I^ N ) T,7=i \J 1 - l^M 
K 2 = 2/{kN)Y j ': = i^ Z ^s 2 /N 2 . 

APPENDIX C 

In order to derive the Laplacian in the respective representation of the matrix field Q, its general definition in an 
arbitrary parametrization, 

A Q = ^Y, dk 9 lk V^d l , (63) 

where the matrix g is the metric tensor, being defined by the quadratic form ds 2 — 1/ATrdQ 2 of the representation 

ds 2 = dx T gdx, (64) 

where x is the vector of parameters of the representation. 

For B ^ 0, Q is element of U(2)/{U{1) x U{1)), by enforcing the conditions Q 2 = 1, Q T C = CQ, and Q+ = Q, 

[Q,n] = 0- 



It can be paramterized as Q — 

where 9 S [0, it] and \ G [0, 2n]. 
Thus, 



cos 9 e lx sin 
e~ lx sin # — cos ( 



ds 2 = d9 2 + sin 2 9dx • (65) 



18 



and g = ( J . ^2 a )• Thus, with Eq. follows: 



A Q = d XD (1 - x 2 D )d Xo + —^rdx 2 , (66) 

1 - A D 

where Xd — cos(0). 

Note that the autocorrelation function depends on the energy difference uj through the coupling TrKQ = 2 * 2 Ad, 
so that only that part of the Laplacian which does not commute with TrKQ, 

A% = d XD (l-\ 2 D )d Xo . (67) 

enters in the frequency dependence of the autocorrelation function of spectral determinants, Eq. ( |3^ ). Since 
U(2)/(U(1) x U(l)) = S2, the two sphere, this is equivalent to the treatment of spherically symmetric potentials, and 
the Laplacian can be identified with the square of the angular momentum, — Aq = L 2 , and L z — id x /(l — Xjj) does 
commute with the Hamiltonian, 

H = -l/(2m)L 2 + ia^^z, (68) 

since z = cos9o does commute with L z . Therefore, u> ^ does not break the azimuthal symmetry of rotations around 
the z-axis, n z . 

For B — 0, Q is element of the symplectic symmetric space, Sp(2) / (Sp(l) x Sp(l)), by enforcing the conditions 
Q 2 = 1, Q T C = CQ, and Q+ = Q. 



One obtains: 



with A = { H * } where I a I 2 + I b I 2 +c 2 = 1 



6* —a" 

A matrix Q with the above symmetries can be represented as, 

Q = U- l QlU, (70) 

with 

U = V C U D , (71) 

where 

U D = V^T^V D , (72) 

where 



cos^c tsmtic T 2 
ism.6cT2 — cos 9 c 



(73) 



and 



and 



o _ ; cos6»d/2 ism9 D /2 . _ 
±D ~ I iwa.0 D /2 cosOd/2 > ' ^ 4J 



and Tj, i = 1, 2, 3 are the Pauli matrices. Such a representation was first given by Altland, Iida and Efetov ]6l[ | to study 
the crossover between the spectral statistics of Gaussian distributed random matrices as the time reversal symmetry 
is broken, within the supersymmetric nonlinear sigma model. Here, in order to study the ASD, we need to consider 
only the compact block of the representation given there. 
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We find that ds 2 = TrdQ 2 /A = d6 2 c + cos 2 c d0 2 D + sin 2 9 C (j) 2 c + cos 2 9 C sin 2 9 D d(j) 2 D and thereby with Eq. (|| ), 
the part of the Laplace operator which does not commute with TrXQ = AXcXd is given by Eq. (32), 

= d Xc (1 - X 2 c )d\ c + 2 l —^8 Xc 

+ l-d XD (l-X 2 D )d Xo) (76) 
A c 

where Aj = cos 9i,i — C,D. 

For moderately strong spin- orbit scattering 1/rgo > Ac, in the functional integral representation of the spectral 
determinants by Grassman vectors the spin degree of freedom is introduced, a — 2 and the matrix C is, due to the 
time reversal of the spinor, substituted by ia>xT\ p2[ |. The spin-orbit scattering reduces the Q matrix to unity in spin 
space. Thus, the matrix C has effectively the form t\. The condition Q T C = CQ leads therefore to a new symmetry 
class, when the spin symmetry is broken but the time reversal symmetry remains intact. Then, Q are 4x4- matrices 
on the orthogonal symmetric space 0(4)/(0(2) x 0(2)) p8[ . 

A matrix Q with the above symmetries can be represented as, 



with 



where 



Q = V-'Q V, (77) 



Q^=( C °1 Sm %), (78) 
1 sm 9 — cos ( 



71 <72 

h 01 



(79) 



with 9i € [0, 7r],« = 1, 2, and 



where 



Vi = exp(ixtT 3 ), (81) 



with Xi € [0,27r],i = 1,2. 

Thus, we find ds 2 = TrQ 2 /4 = Yli=i 2 + ^X T 5xXj where 



_ , sin 2 9i + sin 2 9 2 - sin 2 9i + sin 2 9 2 1 , 
5 * - I - sin 2 9 1 + sin 2 2 sin 2 9 1 + sin 2 2 ' ' 



Thus, the part of the Laplace operator which does not commute with TrXQ = 4AiA 2 is given by Eq. (36), 



A Q= E ^(1-A 2 )^, (83) 
1=1,2 



where Xi = cosf?,, i = 1, 2. 
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